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1. Introduction 

Dyson-Schwinger equations (DSEs) are a nonperturbative continuum approach to 
quantum chromodynamics (QCD). They provide a means to study properties of 
the G reen fu nctions of QCD^I as well as hadrons as bound states of quarks and 
gluons!^™ Hadrons are studied in this framework of infinitely many coupled in- 
tegral equations with the help of a symmetry preserving truncation scheme. In 
the case of mesons discussed here, one solves the Bethe-Salpeter equation (BSE) 
for a quark- antiquark pair. For cal culations of baryon properties, e. g. electromag- 
netic, we ak, and p ionic form factor one uses a covariant set of Faddeev 
equations EE2II3I 

The Dyson-Schwinger equation framework has numerous features, amongst 
them: first, it is a Poincare-covariant framework and thus ideally suited for the 
study of hadron observables such as, e.g., electromagnetic form factors. Secondly, 
symmetries are represented by Ward-Takahashi or Slavnov- Taylor identities, which 
are then built into the scheme used to truncate the infinite tower of coupled integral 
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equations. If a truncation scheme respects such an identity at every step, then one 
can i) prove exact (model-independent) results and ii) use sophisticated models to 
calculate physical quantities which illustrate these results and automatically reflect 
the properties of the corresponding symmetry. 

One such truncation is the so-called rainbow-ladder truncation, which has been 
used extensively and successfully to study meson ground states for more than a 
decade (see, e.g., Refs. El El El El EJ- Despite this success it has become 
obvious that certain states and phenomena, such as axial-vector mesons, exotic 
mesons, a or heavy-light systems are not well-described in the most sophisticated 
calculations availabl e to date. As a consequence, efforts are being made to go be- 
yond this truncation P^W^l^ l ^ lbut these efforts are considerable, forcing present 
sophisticated calculations to remain at an exploratory stage. 



2. QCD Gap and Bethe-Salpeter Equations 

The homogeneous BSE is b 

,A 

[T(p;P)] tu = [x(q;P)] sr K%(p,q;P), (1) 

Jq 

where p is the relative and P the total momentum of the constituents, r,...,u 
represent color, Dirac and flavor indices, 

X (q;P) = S(q + )r(q;P)S(q-), (2) 

q± = qzkP/2, and J represents a Poincare invariant rcgularization of the integral, 
with A the regularization mass-scale. In Eq. Q, S is the renormalized dressed-quark 
propagator and K is the fully amputated dressed-quark-antiquark scattering kernel; 
for details, see Refs. El El 

The dressed-quark propagator appearing in the BSE's kernel is determined by 
the renormalized gap equation 

Sip)- 1 = Z 2 ( l7 .p + m bm ) + S(p), (3) 
S(p) =Zij 9 2 D^(p - q)^- 7fl S(q)T a v (q >P ), (4) 

where D^ v is the dressed gluon propagator, T u (q,p) is the dressed quark-gluon 
vertex, and m bm is the A-dependent current-quark bare mass. The quark-gluon- 
vertex and quark wave function renormalization constants, Zx^^ 2 , A 2 ), depend on 
the gauge parameter, the renormalization point, £, and the regularization mass- 
scale. The leptonic decay constant of a pseudoscalar meson is calculated from the 

a We use the term "exotic" to characterize mesons with quantum numbers that a system composed 
of a constituent-quark and constituent antiquark cannot have. 

b We employ a Euclidean metric, with: {7 M ,7„} = 2<5 M „; 7^ = 7^; 75 = —71727374; 
tr757 M 7„7p7a- = — 4e M „ p(T ; and a ■ b = ^ i=1 Oi&j. For a timelike vector P M , P 2 < 0. 
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solution of Eq. Q via 

frg 5 ij P u = Z 2 tr J -r^ xi> S (q; P) , (5) 
which is gauge invariant, and cutoff and renormalisation-point independent. 



2.1. Rainbow-ladder truncation 

The first step in the symmetry-preserving truncation scheme described in Refs. 1191 
I21ll24ll2"5"l is the rainbow approximation to the gap equation combined with a ladder 
truncation in the BSE. The interaction kernels of Eqs. and © then take the 
form 



K%(p,q;P) = -4ira((f)D%?(Q) 









7„y 


ts 





(6) 



and 



E(p) = | 47ra(Q 2 )I^"(Q)y7 /1 S(g)y7, 



(7) 



where Q = p — q, D*™ C (Q) is the free gauge boson propagator and a(Q 2 ) is an 
effective running coupling. The ultraviolet behavior of this coupling can be taken 
from perturbative QCD, while in the infrared one makes an Ansatz with sufficient 
enhancement to correctly reproduce the phenomenology of dynamical symmetry 
breaking, i.e., enhancement of the quark mass function on a domain p 2 < 1 GcV 2 
and correspondingly a correct value of the chiral condensate. Such an Ansatz i J15 | 16 | 



47ra(s) 



Ait 2 



Dse~ s,u 



In 



3tt 2 7« 



s/A 2 



QCD 



(8) 



with T{s) = [1 - exp(-s/[4m?])]/s, m t = 0.5 GeV, ln(r + 1) = 2, 7m = 12/25 
and Aqcd = ^-jfg = 0.234 GeV. The free parameters of this model are the range w 
and the strength D. One feature of the model is a specific parameter dependence of 
ground-state meson masses and propert ies, namely that they remain constant over 
a range of u>, if oj D = const, is satisfied!^ This results in a one-parameter model, 
where the free parameter uj is varied in the interval [0.3, 0.5] GeV. In a calculation, u> 
and D as well as the current-quark masses are fixed to pion mass and decay constant 
as well as the chiral condensate. That being done, further results are predictions. 



c We use Landau gauge in all calculations. 
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Fig. 1. Pseudoscalar and scalar meson masses for ground and excited states as functions of the 
current quark mass scaled to 1 GeV by one-loop evolution. The three vertical dashed-dotted lines 
(the leftmost is very close to the M-axis) indicate the values for the u/d, s, and c quark masses. 

2.2. Chiral symmetry 

One expression of the chiral properties of QCD is the axial-vector Ward-Takahashi 
identity 

P^(p;P) = 5- 1 (p + ) l7 5y +* 7 5yS- 1 (p_) - 2i m(C) I^(p; P) , (9) 

which is written here for two quark flavors, each with the same current-quark mass: 
{t 1 : i = 1,2,3} are flavor Pauli matrices. T J 5 Jk;P) and T 3 5 (k; P) are the axial- 
vector and pseudoscalar vertices (for details, see Ref. IT5|). Equation © is satisfied 
by relating the kernels of the Bethe-Salpeter and gap equations (JTJ and (J3J, e.g., 
© and |0. A direct consequence of this identity is the relation 

/p S 74 s = 2m(C)p PS (C), (10) 
which relates the pseudoscalar meson mass mps and decay constant fps to the 
current quark mass m(£) and the residue of the pseudoscal ar v ertex at the pion pole 
Pps(C) a t the renormalization point £. It has been showrP^ that the implications 
of Eq. (|10fl in the chiral limit are different for the pion ground and excited states; 
namely, in the presence of dynamical chiral symmetry breaking the mass of the 
ground state vanishes, whereas for the excited states it is the leptonic decay constant 
that vanishes instead. 

3. Mesons 

3.1. Meson ground states 

Ground state mesons and their properties have been studied in rainbow-ladder 
truncation using different Ansatze of various sophistication for the effective run- 
ning couplinj l'^1^'^^01 over a range of quark masses including the heavy-quark 
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Fig. 2. Pseudoscalar and scalar cc meson masses as functions of the model parameter lo. The three 
dotted lines correspond to the experimental values for the lower three states. The r} c ground-state 
mass was fitted to the experimental value to fix the value of the charm-quark mass. 

The Ansatz of Ref. 1161 Eq. ©, has been used success- 
fully to calculate a large number of pseudoscalar and vector meson properties. We 
also used this Ansatz, since it has the correct ultraviolet behavior and thus yields 
reliable results not only for spectroscopy but also for dynamical observables like 
form factors. 

3.2. Radial meson excitations 

It was natural to study radial excitations of pseudo scal ar mesons first. An estimate 
of the excited pion mass and leptonic decay constanlE^and the structure of excited- 
state Bethe-Salpeter am plitu des, which shows similar characteristics to a quantum 
mechanical wave functional were followed by more detailed studies of pseudoscalar 
meson radial excitations ^3 and their electromagnetic properties!^ The present 
work includes a study of scalar mesons and their first radial excitations. 

3.3. Masses 

Figure ^ shows the ground and first radially excited states of pseudoscalar and 
scalar mesons as functions of the current-quark mass. Figs. n an d|3 are generated 
from results for ui = 0.38 GeV. In contrast to the ground states, the masses and 
properties of radial excitations do depend on the value of u>, even if uj D = const. 
Since r = 1/uj corresponds to a range of the infrared part of the interaction, this 
means that radial meson excitations provide a means to study the long-range part 
of the strong interaction. 

To obtain results independent of the choice for u), one can make use of ratios 
of calculated observables, which remain constant over a range of to, to estimate 
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Fig. 3. /ps for the ground- and —fps for the first radially excited-state pseudoscalar as functions 
of the current-quark mass scaled to 1 GeV by one-loop evolution. The three vertical dashed-dotted 
lines correspond to the u/d, s, and c quark masses. Estimated numerical errors are below 5%. 



properties of states on the basis of experimental values of other states that are 
known. An example for such an estimate is that of the mass of the i4f(1460), the 
K radial excitation, where the ratio Mi S r ex /M Wex is calculated to be 1.167; using 
the experimental number M„ ex = 1.3 ± 0.1 GeV this yields Mk S7 . — 1-52 ± 0.12 
Q e y_yy| Meanwhile, we have performed the analogous calculation for the leptonic 
decay constants of th ese sta tes and found fK ex I j-*^ — 10 in agreement with an 
estimate via sum rules l^iilffl 

Figure [21 illustrates the same procedure to estimate the mass of the first scalar 
radial cc excitation Xco(2P). The ratio of the calculated masses for the x c o(2P) to 
the ij c (2S) is 1.066. Via the experimental value of 3.64 GeV (all experimental data 
are taken from Ref. I42J) for the mass of the r) c {2S) we predict the Yrn(2P ) mass to 
be 3.88 GeV. This compares well to quark-model predictionspSEiESEEl which lie 
somewhat below the corresponding estimates from lattice QCD.EZGB1 We note that 
the x c o(2P) has not yet been observed experimentally. 



3.4. Decay constants 

For the pseudoscalar ground and first radially excited states we plot /ps and — /ps, 
respectively, as functions of the current-quark mass in Fig. While /ps is not 
necessarily directly accessible experimentally for QQ systems, where Q is a heavy- 
quark, it is always a well-defined axial- vector moment of the meson's Bethe-Salpeter 
amplitude. Its evolution with current-quark mass is therefore a useful tool with 
which to probe QCD and models thereof. One can see that both curves have a 
maximum at about the c-quark mass; for higher quark masses, the size of /ps 
decreases for both ground and excited states. It is remarkable that this "t urning 
point" occurs around the same quark mass as for the heavy-light case .Isi^lsSI At the 
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Table 1. Values and ratios of /pg for ground- and excited-state pseudoscalar mesons for 
u/d, s, and c quark masses. Current-quark masses are given at the renormalization point 
£ and scaled to 1 GeV by one-loop evolution. All values in GeV. 



Quark 


m(C = 19 GeV) 


m(l GeV) 


/gi 


/ex 


/ex/ /gr 


u/d 


0.0037 


0.00545 


0.092 


0.0015 


0.016 


s 


0.0835 


0.125 


0.13 


0.023 


0.18 


c 


0.905 


1.357 


0.23 


0.19 


0.83 



values for the u/d-, s-, and c-quark masses we extract the values for /ps for the 
ground and ~/ps for excited states as well as their ratio. The results are summarized 
in Tabled We note here that for the s quark our pseudoscalar ground state does 
not correspond to an actual meson, since it consists merely of a n s s component. 
However, a radial ss excitation can be identified with the 77(1475) -ESI Furthermore, 
the ground-state ss properties can also be studied on the lattice, where recent effo rts 
have begun to study r alios of ground- and excited-state leptonic decay constants.^ 

4. Conclusions and Outlook 

We have extended previous studies of radial meson excitations by studying scalar 
excitations and quark masses up to the charm quark. While in the model we used 
ground-state meson properties do not depend on variations of the model parameter 
w, the specific parameter dependence of radial excitation properties allows investi- 
gations of the long-range part of the strong interaction between quarks. Without 
fixing id to a particular value, we used ratios of properties of different excited states 
to make estimates for, e.g., the leptonic decay constant of the i^(1460) and the 
mass of the Xco(2-P)- This is made an efficacious procedure by the fact that ratios 
of excited-state properties remain constant over the domain of ui under investigation 
to a very good level of approximation. For both the ground and excited equal-mass 
pseudoscalar states we have calculated /ps and observed a rise to a maximal size 
around the charm-quark mass and a decrease for higher quark masses. Further ef- 
forts in this direction will include studies of the radial excitations of vector mesons. 
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